Backscattering of Dirac fermions in HgTe quantum wells with a finite gap 
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The density-dependent mobility of n-type HgTe quantum wells with inverted band ordering has 
been studied both experimentally and theoretically. While semiconductor heterostructures with a 
parabolic dispersion exhibit an increase in mobility with carrier density, high quality HgTe quantum 
wells exhibit a distinct mobility maximum. We show that this mobility anomaly is due to backscat- 
tering of Dirac fermions from random fluctuations of the band gap (Dirac mass). Our findings open 
new avenues for the study of Dirac fermion transport with finite and random mass, which so far has 
been hard to access. 

PACS numbers: 



Introduction- Topological insulators, such as HgTe 
quantum wells (QWs) (JQ, Bii_:,Sb^ [1,11, BisSeall, 
and Bi2Te3 0, 0| provide the most recent examples of 
the realization of Dirac fermions in condensed matter 
physics. Unlike graphene [§], where two valleys of Dirac 
fermions exist, in these new materials Dirac fermions ap- 
pear only at a single point in the Brillouin zone. The 
absence of valley scattering and the possibility of un- 
conventional types of disorder, such as a random Dirac 
mass 10], make electron transport studies on these ma- 
terials particularly interesting. However, the low car- 
rier mobility of most topological insulators still poses 
a serious obstacle for transport investigations of spe- 
cific scattering mechanisms. The notable exception are 
MBE-grown HgTe QWs where transport measurements 
have been used to detect the quantum spin Hall state [2] 
and, more recently, to realize massless single-valley two- 
dimensional Dirac fermions |lli] . 

In this paper, we study both experimentally and theo- 
retically a new manifestation of Dirac fermion transport 
in n-type HgTe QWs, i.e., the occurrence of a maximum 
in the mobility as a function of carrier density. The origin 
of the maximum is the competition of two disorder scat- 
tering mechanisms, viz. scattering by charged impurities 
and by QW width fluctuations which induce a fluctuat- 
ing band gap, or, equivalently, fluctuating Dirac mass. 
As in other semiconductor heterostructures [l^l , in HgTe 
QWs the screening of ionized impurities by the carriers 
results, initially, in a monotonic increase of the mobil- 
ity with increasing carrier density. However, while the 
impurity scattering becomes weaker with increasing car- 
rier density, the other scattering mechanism - well-width 
fluctuations - becomes increasingly important, leading to 
a reduction of the carrier mobility. Dirac mass disorder 
generates scattering between states of opposite momenta, 
also called backscattering. Thus, the observed mobility 
peak is a clear manifestation of Dirac fermion backscat- 
tering in HgTe QWs. 

Backscattering of Dirac fermions is most pronounced 
in gapped systems 13|. As a consequence, the mobil- 



ity of graphene does not show a maximum, but rather a 
saturation at high carrier densities (see, e.g. Refs. [qI-ITsI- 
ITsjl). which has been attributed to charged impurity scat- 
tering. Although some theoretical studies indicate that 
mass disorder may play a role in the vicinit y of the neu- 
trality point of graphene (see, e.g. Refs. ]l6l . [l7|). its 
experimental identification has remained problematic. 

Experiment- Transport experiments have been per- 
formed on modulation doped HgTe/Hgo.aCdo.yTe QW 
structures fabricated by molecular beam epitaxy on 
lattice-matched (Cd,Zn)Te substrates. The samples have 
been patterned into Hall bar devices with dimensions of 
(600 X 200) /im^ using a low temperature optical lithogra- 
phy process, and covered by a 5/100 nm Ti/Au gate elec- 
trode which is deposited onto a 110 nm thick Si3N4/Si02 
multilayer gate insulator. Ohmic contacts are provided 
by thermal indium bonding. A micrograph of the struc- 
ture is shown in the inset of Fig. [TJ The samples have 
nominal QW widths, d, ranging from 5.0 to 12.0 nm, 
thus covering both the normal (d < 6.3 nm) and inverted 
{d > 6.3 nm) band structure regimes j^jlUl]- The relevant 
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FIG. 1: Hall and magnetoresistance for sample #6 at zero 
gate voltage. The measurement yields an electron mobility of 
1.1 X lO'' cm"^/Vs and a carrier density of 4.3 x 10^^ cm~^. 
Inset: Micrograph of the Hallbar structure. 
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parameters of all samples are summarized in Table ID We 
have performed standard Hall and Shubnikov-de Haas 
measurements on these samples in magnetic fields up to 
B = 5 T, at a temperature 4.2 K. As an example, Hall 
and magnetoresistance data for sample #6, which has 
the peak highest mobility are shown in Fig. [T] The car- 
rier densities and mobility of ungated samples {Vg = 0) 
are in the range of 3.0 x 10^^ cm~^ to 5.5 x 10^^ cm~^ 
and several 10^ cm^/Vs, respectively. 

The density dependent carrier mobilities are obtained 
from the gate- voltage dependence of the longitudinal re- 
sistance at zero magnetic field, while the dependence of 
the carrier concentration on gate-voltage was deduced 
from the Hall voltage measured at a fixed magnetic field 
of 300 mT. The thick solid lines in Fig. [2] show the ex- 
perimentally observed mobility versus carrier density for 
the various HgTe QWs. In high quality samples with in- 
verted subband structure ordering (#3, #4, and #6) the 
mobility exhibits a distinct maximum for carrier densi- 
ties in the range of 3 to 6 x 10^^ cm~^. For samples with 
a slightly lower mobility (#1, #2, and #5) a saturation 
of the mobility with carrier density is observed. 

Model - In order to explain the unusual dependence of 
the mobility on carrier density we build a model for the 
conductivity in HgTe QWs which is based on the four- 
band Dirac model of Refs. [l|, The effective Dirac 
Hamiltonian has the following form: 

H = s.a-iAk + Mrz) + Dk^ -~E^+Vr + 6M„s,a,, (1) 



where Mg. 



M{d) + Bk?. In Eq. ^ the linear term, 
which is proportional to k = — iVr and the constant 
A originate from the hybridization of the first electron 
(El) and heavy-hole (HHl) subbands in the quantum 
well. These two subbands are represented by the pseu- 
dospin (T, whose components Ux, (Jy and Uz are 2x2 
Pauli matrices. (The real spin degree of freedom is rep- 
resented by the Pauli matrix Sz-) The term proportional 
to the effective Dirac mass reflects the average band 
gap, = |M((i)|, which is determined by the nom- 

inal thickness d of the QW. The k-dependent part of the 
mass (the Bk^ -term) and the parabolic background DkP' 
take further details of the band dispersion in HgTe QWs 
into account [H |3] • Note that the mass term violates the 
pseudo-time reversal symmetry k — )• — k and a ^ — cr of 
Hamiltonian ([1]), which manifests itself in a dependence 
of the conductivity on Mj^. Equation ([1]) also takes the 
two most relevant types of disorder into account: a ran- 
dom potential due to charged impurities, Vr, and spatial 
fluctuations of the Dirac mass, i5Mr, which are related to 
deviations of the QW thickness from its average value d. 
We assume the n-type transport regime, with the Fermi 
level in the conduction band, under the weak scat- 
tering condition, kpV^r 3> 1, where t is the transport 
relaxation time and and A:^ = \/27m are the Fermi 
velocity and wave-vector, respectively. The conductivity 
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FIG. 2: Mobility ^ versus carrier density n for six HgTe QWs 
(#1 ~ #6). The experimental /^(n) dependence (thick lines), 
including the appearance of the maximum, agrees well with 
our model (thin red lines) that takes into account the compe- 
tition of charged-impurity disorder and QW thickness fluctua- 
tions. Inset: Band gap M versus QW thickness d [sample #1 
has a positive band gap, sample #2 has approximately zero 
gap, while the rest ^^3 — #6 are in the inverted regime]. The 
pronounced maximum is observed in high quality inverted 
samples ^^3, #4 and #6 where the impurity concentration is 
sufficiently low (see also Table |I]). 



can then be obtained from the Kubo formula: 
/■ d^k 

"2^ 



rTr 



(2) 



where G^'e disorder-averaged retarded and advanced 
Green's functions and Vx, Vx are current vertices (the 
tilde indicates vertex renormalization by disorder in the 
ladder approximation [13]). The resulting conductivity 
axx is proportional to the density of states per spin at 
the Fermi level, = kp /inhvp, and the transport time 
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T = 1/ Jq^ {I — cos9)w{9)d9, where w{0) is the scatter- 
ing rate at angle 9. Below, we obtain an expression for 
w{9) from the electron self-energy in the self-consistent 
Born approximation (SCBA), calculate r and the carrier- 
density dependent mobility ^(n) = CT2:2;(n)/ne. 

For this purpose, we assume that the potential and 
mass disorder are uncorrelated and completely charac- 
terized by the two-point correlation functions (T4^r') = 
C^r' and {SMr6Mr') = C^^' • In k space this leads to the 
Dyson equation Gj^. = Gn 



^ok^k^k fo"" the disorder- 



averaged Green's function where the self-energy T,^ 
is given by the standard SCBA expression 



^ok — - — 
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The unperturbed Green's function Gq^. describes 
a conduction-band electron with dispersion ^k = 
^^^2^2 + + Dk^ - and chirality s^a- ■ e^, = 1 



(/ is the unit matrix). The solution [IE 
function = ^(/ -I- SzCr ■ ek)/(e - fk 

the finite elastic life-time t^i — 1/ J^^ w{9)d9, where the 
scattering rate w{9) at angle 9 is given by 



for the Green's 
-I- 27^) contains 



N 

MO) = -f 



(4 



I sine/2| 



s/2|ei } cof? - 



^2k J-, sin 



e/2\^l- + C2fc^|sm6'/2| ' ^^^^ ~ 



(5) 



Here e_L = Mkp/ ^ A'^k^ + M"^^ is the out-of-plane com- 
ponent of the unit vector Ck at |k| = fc^ [see Eq. (j4|)]. 

Analyzing Eq. ([SJ, one notes that the first term van- 
ishes at ^? = TT. This is the well-known absence of 
backscattering in the limit of massless Dirac electrons [IH 
(this behavior is plotted as the red curve in Fig. [3K)- In 
gapless Dirac materials, the pseudospin cr points along k, 
and, therefore, states with k and — k are orthogonal to 
each other, i.e. unavailable for scattering. The most es- 
sential distinction between HgTe quantum wells and the 
zero gap case is the second term in Eq. ([5|), which actually 
has a maximum at the backscattering angle 9 = ir. This 
term originates from the finite Dirac mass Mj^^SzCz and 
its spatial fluctuation dM^SzCTz- Both of these result in 
an out-of-plane pseudospin component <7z , hence the op- 
posite k states are no longer orthogonal and large-angle 
scattering now becomes possible (cf. the black curve in 
Fig. [3^). Figure [3b shows that the backscattering is en- 
hanced with increasing carrier density n, which accounts 
for the non-monotonic behavior of fi{n) in inverted quan- 
tum wells observed in Fig. [2l 

To proceed further, we make specific assumptions for 
the correlation functions: 
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FIG. 3: (a) Polar plot of the angular dependence of the 
normalized scattering rate wiO) for massless and massive 
(M = -15 meV and C,'^'' = 1.5 x 10"^) Dirac fermions [see 
Eqs. ((5]) and Massive carriers show strong backscattering 
at angles 7r/2 < 6 < 37r/2. (b) Polar plot of the normalized 
function w{9) for different carrier densities n in the inverted 
regime (M = -15 meV and = 1.5 x 10"^). The large- 
angle scattering is enhanced with increasing n, which accounts 
for the non-monotonic mobility ^i(n) in Fig. (2] 



where is the usual correlation function of screened 
Coulomb impurities [isl [l4| with the concentration 
and the average dielectric constant, e w 15, of the 
HgTe/CdTe QW. Ck is normalized such that C*^ is a 
small dimensionless parameter, which guarantees that 
SM is small compared with the leading linear term 
Akp in Hamiltonian ([T]). Since 5M is caused by fluc- 
tuations of the QW thickness, C*^ is independent of car- 
rier density. As shown below, this approximation yields 
an excellent agreement with the measurements. Using 
Eqs. ^ and ([5]) we find the mobility = (T^x/ne: 



where 



4e//i 



M [3+cl{n)]A- 
h?v^ (n) 



-,(7) 



ej. (n) = {M + 2tt Bn) / ^/2TTA^n + (M + 2TrBn)^ , (8) 
CY = d9sm'^9[l + ehv^{n)e-^sm9/2]^^ , (9) 

/■TT 

CY / d9 sin^ 9 /2[l + ehv^{n)e-^sm9/2]^^ .{10) 
Jo 

Using Eq. ^ we can quantitatively reproduce all ex- 
perimental curves for in Fig. [5] using the disorder 
parameters n^, (^^^ and band gap M indicated in Table H] 
The values of M obtained from this fit agree well with 
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TABLE I: Sample and disorder parameters (the other param- 
eters of the model are fixed: e = f5, A — 0.38 eV-nm, B = 1.2 
eV-nm^ and D = 0.85 eV nm^). 

#f#2#3#4#5 #6 

nominal well width (nm) 5.7 6.3 7 7 7.5 12 

max. mobility (fO^ cmVVs) f.32 2.96 3.71 4.76 3.33 12.27 

M (meV) 10 -12 -12 -15 -24 

ni(xlO^°cm-^) 8 3.41 2.95 2.55 4.56 1.09 

C*^(xlO-^) 1.65 0.5 1.65 0.91 0.43 0.33 



those obtained from band structure calculations and the 
analysis of the experimental SdH oscillations. 

Discussion.- Using the theoretical model presented in 
the previous section, we can now explain the observed 
non-monotonic dependence of /i(ri) as follows. The ini- 
tial increase in results from the impurity screen- 
ing: it reflects the density-dependence of the Fermi ve- 
locity (n) which enters the screened impurity poten- 
tial via the DOS in Eq. (H]) [see also Eqs. © and pi!]) ]. 
Similar behavior was found in conventional doped het- 
erostructures fl^ . However, here Vp{n) is specific to the 
massive Dirac Hamiltonian ([T|). Furthermore, for the in- 
verted quantum wells #3 — #6, the mobility is addition- 
ally enhanced due to the reduction of the total Dirac gap 
— |Af((i)| + 2'KBn at low carrier densities n [cf. Eq. (|5])]. 
This leads to a more rapid initial increase in com- 
pared to sample ^1, which has a normal band structure, 
and the zero-gap sample ^2. At higher carrier densities 
the mobility starts to decrease for all the inverted sam- 
ples, most pronouncedly so for the high-quality samples 
#3, #4 and #6. Since the estimated impurity concentra- 
tion is lowest for these samples (see Table we attribute 
this decrease to the fluctuations of the QW width (Dirac 
mass) , accounted for by the term cx nC*^ in Eq. ([7]) . Intu- 
itively, the reduction of the mobility can be explained by 
the fact that the rate of scattering off the well-width fluc- 
tuations grows proportionally to the carrier DOS (n) 
because more states become available for backscattering 
as the Fermi surface size increases with n [see Eqs. (O 
and dSI]. 

From the fits we estimate the amplitude of the well 
width fluctuations, Sd, to be of the order of 0.2 — 
0.3 nm, which is obtained by integrating the correla- 
tion function of the thickness fluctuations, {SdrSdr') — 
{6 Mr6 Mr') / M' (d)'^ , over the area. This integral is a mea- 
sure of the typical height 6d times length L of the fluctu- 
ation: L6d ~ A^/C^/\M'{d)\ « 0.7 nm^, where we take 
A = 0.38 eV-nm and = 10"^ from Table U and the 
gap derivative |M'(o?)| w 17 meV-nm^^ from the inset in 
Fig. [2j For a realistic sample L ^ 6d, thus we estimated 
the ratio L/5d 10, which yields 6d ~ 0.26 nm. The re- 
sult is in good agreement with X-ray reflectivity data on 
MBE-grown HgTe QW structures of similar quality [l9| . 

Conclusions.- We have shown both experimentally and 



theoretically that the density-dependent mobility of high 
quality HgTe quantum wells with an inverted band struc- 
ture exhibits a maximum. While the initial increase 
in the mobility is mainly due to scattering at screened 
charged impurities, the decreasing part is associated with 
the band gap fluctuations that generate mass disorder 
for Dirac-like fermions in this material. Our flndings 
thus clearly demonstrate the occurrence of Dirac fermion 
backscattering in finite gap systems. 
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